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Abstract 

Optimum trajectories for the Applications 
Technology Satellite (ATS) -E mission are obtained. 
Analysis, procedure, and results are presented. 

The trajectories are numerically integrated from 
launch to insertion into the final orbit. As a 
result of a much smaller than optimum apogee motor, 
these trajectories, unlike conventional synchro- 
nous orbit trajectories, require non-circular park- 
ing orbits and large amounts of inclination reduc- 
tion before the solid motor burn at apogee. Con- 
straints on parking orbit perigee radius and dura- 
tion are included. Figures are presented describ- 
ing the results . 

Introduction 

The Applications Technology Satellite (ATS)-E 
mission Is a circular synchronous equatorial orbit 
mission. The ATS program has the objective of 
advancing technology In areas which may have appli- 
cation to future spacecraft. The experiments which 
are conducted are spacecraft, communication, and 
science oriented. 

The space craft -or iented experiments on the 
ATS-E provide information on power supply and con- 
trol systems, a gravity -gradient stabilization 
system, resistojet and Ion micropound thrusters, 
and synchronous environment. The scientific- 
oriented experiments gather data on the particle 
(electron and proton) distribution and flux and 
the character of the electric and magnetic fields 
at synchronous altitude. 

The launch vehicle for the ATS-E mission was 
an Atlas -Centaur and the solid apogee motor was a 
part of the spacecraft system. The apogee motor 
total impulse was sized for the early ATS missions 
on the Atlas -Agena launch vehicle, which has less 
payload capability than the Atlas -Centaur . The 
apogee motor, although smaller than optimum for 
the larger vehicle, remained unchanged. 

For an optimally sized apogee motor, a con- 
ventional trajectory to circular synchronous 
equatorial orbit is near optimum. A conventional 
trajectory consists of five consecutive phases 
as shown in Fig. 1. The first phase Is an ascent 
from the launch site to a circular parking orbit. 

To maximize the mass In orbit, a 90° launch azi- 
muth is used, which results in a parking orbit 
inclination equal to the launch site latitude. 

This inclination must be removed during the tra- 
jectory. The second phase Is a coast arc to the 
proximity of the equator. A small portion of the 
required plane change is removed by the second 
burn, the third phase. Much more Importantly, the 
second burn must place the vehicle In a transfer 
orbit whose apogee Is over the equator and equal 
to synchronous altitude. The vehicle coasts to 
apogee in the fourth phase. The fifth phase con- 
sists of a final burn that removes the major por- 
tion of the Inclination and circularizes the orbit. 


Intuitively, It seems reasonable that this conven- 
tional profile is near optimum if the burn and 
coast durations may be varied to maximize the mass 
at the end of each burn. However, If the total 
Impulse of the final burn is fixed at less than 
the optimum value, the conventional trajectory 
must be modified to yield maximum payload to the 
final orbit. In particular, the parking orbit is 
noncircular, the perigee radius of the transfer 
orbit increases, and the second burn removes more 
than a minor part of the inclination. The optimi- 
zation problem is to find the best combination of 
these changes and other less important ones to 
yield maximum payload to circular synchronous equa- 
torial orbit. 

Optimization of the conventional trajectory to 
circular synchronous equatorial orbit has been , . 
treated by several authors. Hoelker and Silber' ' 
present a detailed analysis of the conventional 
problem. Rider' ' considers the problem of chang- 
ing the plane and also the radius of a circular 
orbit. These and other similar studies treat the 
problem as one of changing the plane and radius of 
a circular orbit, ignoring the ascent to the first 
(parking) circular orbit. This is satisfactory 
for the conventional case. However, an unconven- 
tional trajectory is more complex since the parking 
orbit Is in general noncircular. The ascent must 
be Included as part of the optimization problem. 
Therefore, a more sophisticated optimization pro- 
cedure Is required for unconventional trajectories. 
Additionally, the references mentioned above are 
general and consequently are not concerned with 
constraints which may alter the acceptability of a 
given trajectory, such as limitations on coast time 
or the minimum perigee radius of the noncircular 
parking orbit. 

The problem of optimizing trajectories to 
circular synchronous equatorial orbit may be con- 
sidered as a multistage launch vehicle optimization 
Several analyses have been performed to optimize 
multistage launch vehicles, including one by the 
authors of this report^ ^ . For optimizing the un- 
conventional trajectory, the analysis in Ref. 3 
was expanded to three dimensions and also, to In- 
clude a constraint on the parking orbit perigee 
radius. The perigee radius constraint must be 
Included to limit aerodynamic heating on the space- 
craft. 

The Applications Technology Satellite (ATS)-E 
mission on the Atlas -Centaur vehicle requires an 
unconventional trajectory to achieve maximum pay- 
load. The final burn is performed by a solid 
motor which is part of the spacecraft system. That 
motor is significantly smaller than optimum. 

There are spacecraft and launch vehicle constraints 
on the trajectory which must be incorporated into 
the solution. The perigee radius and the parking 
orbit coast duration are limited. The results for 
this mission are presented. 
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Analysis 

Problem Description 

A conventional trajectory to circular syn- 
chronous equatorial orbit launched from the Eastern 
Test Range consists of five phases. They are: 

1. Ascent to parking orbit. 

2. Parking orbit coast. 

3. Second impulse. 

4. Transfer orbit coast. 

5. Third impulse or apogee burn. 

Figure 1 shows the planar characteristics of 
the conventional trajectory. The nonplanar charac- 
teristics are shown in Fig. 2. The vehicle is 
launched at an azimuth of 90° in order to maximize 
the vehicle mass in parking orbit and to minimize 
the inclination of the parking orbit. The circular 
parking orbit altitude is as low as aerodynamic 
heating constraints will allow, usually about 165 
to 185 kilometers. The parking orbit coast time 
is usually about fifteen minutes - the time re- 
quired to coast from orbit insertion to the first 
equator crossing. The third phase places the 
vehicle in a transfer orbit whose apogee and peri- 
gee are over the equator. The apogee altitude is 
about equal to the required altitude for a circular 
synchronous orbit. The transfer orbit coast time 
is about five and one-half hours. The third im- 
pulse, the apogee burn, occurs at apogee of the 
transfer orbit. Apogee is designed to occur at the 
equator and at the proper altitude for injection 
into the final orbit. A small part of the inclina- 
tion is removed by the second impulse with the re- 
mainder being removed by the apogee burn. In this 
conventional method, the final conditions at the 
end of each burn are known and the mass can easily 
be maximized progressively phase by phase if the 
second and third impulse sizes are unspecified. 

Now suppose that the total impulse of the 
third burn is fixed. Then the transfer orbit must 
be constructed such that the AV available from the 
third Impulse Is exactly that required to place 
the vehicle In circular synchronous equatorial 
orbit. If the AV available from a fixed third 
total impulse Is less than that required to cir- 
cularize and equatorialize the orbit for the mass 
available from a conventional ascent and second 
Impulse, then the trajectory to transfer orbit 
insertion must be altered to reduce the AV required 
of the third impulse. This can be done by reducing 
the required plane' change and the AV required for 
circularization. For reasons described at length 
in the Results and Discussion section, the uncon- 
ventional trajectory needed to reduce the AV re- 
quired of the apogee motor varies in many respects 
from the conventional profile. The most dramatic 
changes are a noncircular parking orbit, nontrivial 
inclination reduction by the second impulse, and a 
significantly nonequatorial latitude for the second 
impulse. As is desired, the changes result in 
lowering the AV required of the fixed apogee motor. 
However, in this unconventional profile, the final 
conditions required at the end of the ascent and 
second Impulse are unknown. They might be deter- 
mined by varying those final conditions parametri- 
cally until the optimum is obtained. However, 
because of the number of variables, this process 
is clumsy and time consuming. 


Calculus of Variations Solution 

A Calculus of Variations formulation was used 
to maximize the payload to circular synchronous 
equatorial orbit without resorting to a parametric 
search. The optimization of the atmospheric por- 
tion of the trajectory Is omitted from the varia- 
tional analysis since the steering is constrained 
by factors other than optimizing performance, such 
as aerodynamic loading and heating limitations. 

The analysis considers the problem from the point 
in the trajectory that the atmosphere can be ne- 
glected to insertion Into the final orbit. In 
addition to optimizing the steering, the durations 
of any unspecified burns and coasts are optimized 
while maintaining the specified perigee radius of 
the parking orbit. The analysis is presented in 
appendix B. It is derived in three dimensional 
rectangular coordinates In a manner similar to 
Ref. 4. The equations for optimum burn and coast 
duration are obtained from an analysis similar to 
that used by the authors in Ref. 3. It is neces- 
sary to extend the analysis to include an inter- 
mediate boundary condition which specifies the 
perigee radius of the parking orbit at the end of 
the ascent. Additionally, the oblate earth model 
must be added to the variational analysis. The 
effect of oblateness Is not negligible in trajec- 
tories to circular synchronous equatorial orbit. 
Trajectories to that orbit are long, minimally 
around six hours. Oblateness Is the major perturb- 
ing force during most of a trajectory. Because of 
the large change In inclination required to perform 
the mission, any perturbation in the Inclination, 
thus Increasing or decreasing the amount of plane 
change required of the propulsion systems, affects 
the final mass and should be considered in the 
analysis. 

The trajectories are numerically integrated 
to incorporate a nonimpulsive vehicle model and to 
Include the effects of oblateness and small thrusts 
over long periods of time which cannot be conven- 
iently treated impulsively. 

The analysis presented in appendix B requires 
the solution of a two point boundary value problem. 
The solution to the two point boundary value prob- 
lem for the circular synchronous orbit problem with 
a fixed apogee burn and parking orbit coast time re- 
quires satisfaction of a minimum of eight final con- 
ditions with an equal number of Initial conditions. 
The number and specific initial and final conditions 
are explained in appendix B. 

Procedure 

A simple Newton Raphson iteration scheme was 
used to solve the two point boundary value problem. 
This scheme was used successfully with as many as 
twelve Iteration variables. For further explana- 
tion of the Iteration scheme, see Ref. 5. 

The partial derivatives required for the iter- 
ation scheme were obtained by integrating the ad- 
joint equations. These were obtained as in Ref. 4. 
Solutions were initially obtained by using a spher- 
ical earth model for the adjoint equations, but it 
was found that including the oblateness terms im- 
proved the convergence properties of the problems. 

In some problems of this type, It was found that 
Including the oblateness terms was necessary to 
obtain convergence. 
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It was difficult to obtain solutions to these 
problems because of the high degree of nonlinearity 
of many of the derivatives as well as the diffi- 
culty of guessing at the initial values of the 
thrust angle in pitch and yaw and their rates* A 
technique was devised to systematically proceed 
from a simple, easily converged problem to the 
final solution. This technique is described at 
length in appendix C. Other techniques, such as 
gradient methods, might avoid some of the diffi- 
culties associated with the Newton Raphson techni- 
que. However, the method described in the appendix 
Is convenient, straightforward, and adequate. 

After obtaining one solution, proceeding to others 
in the region of interest Is not difficult. 

Results and Discussion 
I . La unch Vehicle 

The spacecraft Is launched by an Atlas - 
Centaur , a two -and -a -half stage vehicle. The Atlas 
is propelled by two booster engines and one sus- 
tainer engine. The booster engines are jettisoned 
at a predetermined acceleration level. The sus- 
tainer engine continues to burn (sustainer solo). 

The Centaur insulation panels and then the payload 
fairing are jettisoned In this phase. The sus- 
tainer solo ends at propellant depletion and the 
Atlas stage is jettisoned. After about ten seconds, 
the Centaur engines, burning hydrogen and oxygen, 
Ignite and burn until the desired parking orbit Is 
reached. During the parking orbit, a hydrogen 
peroxide propulsion system Is used to maintain a 
very small acceleration for propellant retention 
and for attitude control. At the end of the park- 
ing orbit, the dentaur engines burn again until 
the proper transfer orbit Is achieved. After en- 
gine shutdown, the Centaur control system acquires 
the proper orientation for the spacecraft burn, 
the Centaur and the spacecraft separate, and the 
spacecraft is spun up for stability. The space- 
craft coasts up to the proper altitude maintaining 
the separation attitude. The spacecraft motor 
bums to place the spacecraft in the final orbit. 
The spacecraft apogee motor has thrust of 22 240 
newtons and an effective specific impulse of 
279.1 seconds. The total Impulse available from 
the motor is 950 900 newt on -seconds, which corre- 
sponds to a propellant load of 347 kilograms. 

II. Trajectory Description 

The trajectory starts with a short vertical 
rise, followed by a rapid pitchover phase in the 
desired azimuth direction. The amount of pitch- 
over determines the amount of lofting during the 
atmospheric portion of the trajectory. The re- 
mainder of the atmospheric phase (which is assumed 
to end at booster stage jettison). Is flown with a 
near -zero angle of attack steering program (des- 
cribed in Ref. 5), to minimize vehicle heating and 
aerodynamic loads . The thrust direction Is con- 
strained to be parallel to the launch azimuth 
plane, which Is established at launch. 

Since the ATS-E spacecraft motor has a fixed 
propellant load, the trajectory must be designed 
such that the AV required at apogee of the transfer 
orbit is exactly that required to place the space- 
craft in the desired final orbit. As mentioned 
earlier, the ATS-E motor Is much smaller than opti- 
mum. The Atlas-Centaur can put more mass in a con- 


ventional transfer orbit than the apogee motor can 
place in circular synchronous equatorial orbit. 
Therefore, an unconventional trajectory is required 
to lover the A V required of the apogee burn. 

An optimum unconventional trajectory was ob- 
tained for the ATS-E mission to circular synchro- 
nous equatorial orbit. The AV required of the 
apogee motor is reduced by decreasing each of the 
two components which together make up the total AV - 
that needed to circularize the orbit and to reduce 
the inclination to zero. The AV for circulariza- 
tion is. reduced by increasing the horizontal velo- 
city at apogee of the transfer orbit without adding 
radial velocity. Any radial velocity would have to 
be removed by the apogee burn. Increasing the 
horizontal velocity at a fixed apogee radius is 
equivalent to raising the perigee radius of the 
transfer orbit - thereby decreasing the ellipticity 
of the transfer orbit. 

The AV required at apogee for reducing the 
inclination to zero is decreased by lowering the 
inclination of the transfer orbit. However, 
raising the velocity at apogee increases the AV 
required for Inclination removal at a fixed trans- 
fer orbit inclination. Therefore, the combination 
of the two methods represents a compromise which 
Is optimized as part of the total problem. 

In order to obtain the modified transfer orbit, 
the trajectory to insertion into that orbit Is 
modified. Most of the inclination reduction is 
performed by the second burn near the equator. 

Only a small part of the Inclination change to 
transfer orbit Insertion is accomplished in the 
ascent to parking orbit. 

The characteristics of the optimum parking 
orbit are changed from the conventional profile to 
increase the perigee radius of the transfer orbit. 
An elliptical rather than circular parking orbit 
is used to raise the altitude of the second burn. 
The perigee radius of the optimum parking orbit 
remains limited by aerodynamic heating considera- 
tion at some acceptable value. Since injection 
into the parking orbit occurs near perigee, the 
vehicle must coast along the ellipse to a higher 
radius. Due to limitation of the coast duration 
for the ATS-E mission, the second burn was required 
to occur near the first equator crossing. (From 
tracking or other considerations, a second (or 
greater) equator crossing could be chosen for the 
second burn, which would increase the parking or- 
bit coast time by a half period (or more)). The 
latitude of the second burn is no longer equatorial 
as in the conventional, case since the optimum posi- 
tion for raising the perigee radius and decreasing 
the inclination is dependent on radius and velocity 
as well as latitude. The parking orbit coast time 
is greater for this unconventional profile since 
the time to the equator is greater for an ellipti- 
cal than for a circular parking orbit and addition- 
ally, the second burn occurs significantly south 
of the equator. Optimum true anomalies are found 
for the beginning and end of the parking and trans- 
fer orbit coasts. In addition, the optimum com- 
bination of the changes just described as charac- 
terizing the unconventional profile is selected. 

The desired final inclination for the ATS-E 
mission is not exactly zero. The perturbations of 
the sun, moon, and oblateness of the earth cause 
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a spacecraft to drift from an exactly equatorial 
orbit. Since zero inclination is not a stable 
condition , a final orbit inclination yielding the 
smallest average inclination over the lifetime of 
the spacecraft is desired. Small final inclina- 
tions with the proper inertial ascending node are 
found to yield acceptable inclination over the 
lifetime of the satellite. The particular combi- 
nations of final orbit inclination and ascending 
node are functions of the positions of the sun and 
moon, which are in turn functions of launch time 
and date. Therefore, data were obtained for pay- 
load to circular synchronous orbits as a function 
of final inclination. Negative inclinations are 
included in the data. This convention indicates 
that the node has been switched approximately 
180° by the apogee burn. 

The Atlas -Centaur has a twenty-five minute 
limitation on parking orbit coast time for the 
mission. Therefore, inclusion of that constraint 
is necessary for realistic determination of vehicle 
capability. However, optimizing the coast time 
provides a more dramatic and obvious demonstration 
of the optimization procedure. Launch azimuth was 
not optimized along with the other trajectory 
parameters. The effect of launch azimuth was 
investigated parametrically to determine its ef- 
fect on separated spacecraft mass. 

Figure 3 presents separated spacecraft mass 
as a function of launch azimuth for final inclina- 
tions of (-)2° and 5.25° for both optimum and 
twenty-five minute parking orbit coast times. 
Separated spacecraft mass is the mass of the 
spacecraft when it is separated from the Centaur 
vehicle. This figure shows that the separated 
spacecraft mass is rather insensitive to launch 
azimuth. Hence., for simplicity, launch azimuth 
is fixed at 90° for the remaining figures. 

Figure 4 shows the separated spacecraft mass 
as a function of final inclination. The separated • 
spacecraft mass decreases as final inclination 
decreases. Figures 5 and 6 show the effect of 
final inclination on the transfer orbit inclina- 
tion and inertial velocity at apogee. As might 
be expected, as the final inclination decreases, 
so does the transfer orbit inclination. 

As might not be expected, the velocity at 
apogee also decreases as final inclination de- 
creases. Figures 7 through 12 show why this 
occurs. Figure 7 shows the latitude of the second 
Centaur engine start as a function of final in- 
clination. Since the second burn is required to 
remove more inclination as final inclination de- 
creases, it is advantageous to move the burn 
nearer the equator for more efficient plane change. 
Figure 8 shows that the longitude of second burn 
start also decreases as final inclination de- 
creases. These trends decrease the parking orbit 
coast arc as final inclination decreases. This 
is reflected in a decrease in the true anomaly at 
second Centaur cut-off, as seen in Fig. 9. Fig- 
ures 10, 11, and 12 also show additional effects 
of moving the second burn nearer the equator. It 
decreases the parking orbit coast time, the alti- 
tude of the second burn, and the apogee altitude 
of the parking orbit. These all occur as a result 
of the decrease in parking orbit coast arc. These 
figures show why the apogee velocity is decreasing 
as final inclination decreases. The perigee 


radius of the transfer orbit decreases as the alti- 
tude of the second burn decreases. The apogee al- 
titude of the transfer orbit is almost constant at 
synchronous altitude, hence as perigee decreases, 
so does apogee velocity. 

Figures 5 and 6 also indicate that more AV is 
required of the apogee motor as final inclination 
decreases. It can be seen that both the plane 
change and circularization AV are increasing. How- 
ever, Fig. 4 shows that the ignition mass of the 
fixed solid motor is decreasing, which increases 
the AV capability of the apogee motor. 

Figure 13 shows the percentage of the Centaur 
propellant used in the first burn. The figure 
shows that as the final inclination increases, the 
first burn duration increases as the apogee alti- 
tude increases (Fig. 12). 

The final longitude as a function of final 
inclination is shown in Fig. 14. It shows that 
longitude decreases as final inclination increases. 
The satellite remains at the longitude indicated 
only when the inclination is zero, the orbit cir- 
cular, and the altitude synchronous. For other 
inclinations, the position (latitude and longitude) 
of the satellite subpoint describes a figure eight 
on the surface of the rotating earth. The longi- 
tudes indicated in Fig. 14 are injection longitudes, 
not necessarily the longitude at which the equator 
crossing occurs. For small inclinations, the lon- 
gitude does not vary greatly during the period of 
the orbit. 

Now consider the limitation of parking orbit 
coast time. Twenty-five minutes is less than 
optimum for all the final inclinations considered, 
as seen in Fig. 10. The differences in separated 
spacecraft mass are shown in Fig. 4. As seen from 
these figures, as the difference between the opti- 
mum and limited coast times decreases, the loss in 
payload due to coast time limitation decreases also. 

The coast time limitation reduces the advan- 
tage of raising the apogee of the parking orbit as 
final inclination increases. The energy required 
to raise apogee does not yield the payload in- 
creases available with optimum coast time since 
the altitude cannot be acquired as efficiently in 
the shorter coast time. The energy is better 
spent by the second burn to reduce the inclination 
of the transfer orbit. This is reflected in sev- 
eral of the figures. In Fig. 5, the transfer or- 
bit inclination for the coast limited case lies 
well below the optimum case. The lower second 
burn altitude is reflected in the lower velocity 
at apogee of the transfer orbit, as seen in Fig. 6. 
Because the parking orbit characteristics do not 
vary greatly with final inclination, the latitude 
and longitude of the second burn and the true 
anomaly at second Centaur cut-off are nearly con- 
stant. These may be seen in Figs. 7, 8, and 9. 

The conclusions which may be drawn for the per- 
centage of Centaur propellant used in the first 
burn and final longitude (Figs. 13 and 14), are 
similar to those for the optimum coast case. 

Summary of Conclusions 

Analysis and results are presented for tra- 
jectories to circular synchronous equatorial orbit 
where the apogee motor is fixed at a smaller than 
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optimum total impulse. The results for the small 
apogee motor case were obtained for the ATS-E 
mission, which used the Atlas -Centaur launch ve- 
hicle. 

The results show some of the characteristics 
of optimum trajectories for launch vehicle -apogee 
motor combinations where the apogee motor is smal- 
ler than optimum. More important, the results 
demonstrate that optimum trajectories to circular 
synchronous equatorial orbits may be obtained with 
detailed and hence complicated vehicle models for 
unconventional (small apogee motor) trajectory 
profiles. These results may be obtained without 
resorting to exotic mathematical procedures for 
solving the two point boundary value problem. 

These results were obtained with a simple Newton- 
Kaphson iteration scheme. The partial derivatives 
were obtained by integrating the adjoint equations. 
The simple iteration scheme with the integrated 
partial derivatives is able to obtain solutions to 
the highly nonlinear two point boundary value 
problem even when the number of initial and final 
conditions reaches twelve. 

Appendix A 

Symbols 

C first integral of Euler -Lagrange equa- 

tions, kg/sec 


z unit vector pointing at north pole, N.D. 

(3 mass flow rate, kg/sec 

e Jump factor 

T) Lagrange multiplier, kg/sec 

X Lagrange multiplier, kg-sec/m 

li Lagrange multiplier, kg/m 

a Lagrange multiplier, N.D. 

pitch attitude, keg. 
q> yaw attitude, deg. 

Superscripts 

* time derivative 

vector 
unit vector 
f final 

o initial 

Subscripts 


e 

E 


f 

F 

g 

G 


G l> G 2 


h 

J 

m 

N 

P 

r 

r P 

s 

t 

T 

v 

X 


eccentricity, N.D. 

2 / 7 

energy per unit mass, m /sec 

unit thrust direction, N.D. 

functional defined by equation (3) 
kg/sec 

intermediate boundary equation 
gravity acceleration, m/sec^ 


i ) J ) 

£,m,n stage numbers 

f final 

o initial 

d desired 

pk parking orbit 

Operators 


components of oblate gravity accelera- 
tion, m/sec 2 

spherical earth gravity constant, m 3 /sec 2 
angular momentum per unit mass, m /sec 
functio nal to he minimized, kg 
mass, kg 

total number of stages, N.D. 
semi-latus rectum, m 
radius, m 
perigee radius, m 

variational switching function, N.D. 
time, sec 
thrust, N 

velocity, m/sec 
state variable 


dot product 
cross product 
differential 

gradient with respect to x 
partial derivative 

Appendix B 
Derivation of Optimum Control 

As mentioned in the Analysis Section, the 
optimization of a trajectory to a circular synchro- 
nous equatorial orbit may be considered as the pro- 
blem of optimizing a multi-stage launch vehicle to 
a particular final orbit. The optimization problem 
to be considered here begins at booster Jettison, 
which is assumed to be a fixed position and velo- 
city. The sustainer portion of the Atlas continues 
until propellant depletion. The sustainer is jet- 
tisoned and a few seconds later the first Centaur 
burn begins. Its duration is variable and must be 
optimized. The perigee radius of the parking orbit 
is fixed. The duration of the parking orbit 
may or may not be optimized. The parking orbit is 
not a true coast since a small acceleration is 


d( ) 
*( ) 
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maintained for propellant retention. The duration 
of the second Centaur burn must be optimized, fol- 
lowed by an optimum transfer orbit coast (a true 
coast), and a final burn of fixed total impulse. 

The analysis presented in this appendix to solve 
this problem is a special case of the analysis 
derived in Ref. 3, with an additional constraint - 
parking orbit perigee radius. 

The variational problem to be solved is to 
find the steering program and various stage dura- 
tions which maximize the payload capability of a 
multi -stage launch vehicle to a specified final 
orbit. The trajectory must satisfy certain ini- 
tial, final, and intermediate conditions on the 
state variables. The thrust, propellant flow rate, 
and jettison weight for each stage are assumed to 
be constant. The equations of motion and con- 
straints for each stage may be written as 


The optimum thrust direction f is obtained by com- 
bining equations (id) and (5d) and using the 
Ueierstrass E-test. This procedure results in 

f - A (6) 


Integrals of the Motion 

Since F does not explicitly depend on time, 
an integral of the motion is 


C + A • G + IT • v + 1 A - cjPi =0 (7) 

When a spherical gravity model is assumed (i.e., 
G(r) * Gm/r 3 r), three additional integrals of 
the motion exist which are given by 


AXv+qXr= constant 


V - G (r) -4 f = 0 

(la) 

r - v = 0 

(lb) 

m + = 0 

(lc) 

A A 

f ■ f - X = 0 

(Id) 

A 

where f is the unit thrust direction and G 
the oblate ’earth gravity acceleration^, 
which may also be written 

(r) is 

(j(r) = G 1 (r, r ■ z) r + G 2 (r, r * z)z 

: (2) 

(All symbols are defined in appendix A.) Suppose 
that each stage of the vehicle is numbered con- 
secutively starting with the booster. For analy- 
sis purposes a stage change occurs when the thrust 
and/or propellant flow rate changes and/or a mass 
is jettisoned. A Bolza formulation of the varia- 
tional problem is usedA^ , and the functional to 
be minimized is written as in Ref. 3 as 

N f t i 

J = -m- dt 

t i-l 

(3) 

where the functional F i for each stage is 


F ± = A -[v - 0 - ^-f] + 7 • [? - ▼] 


+ a(m + 0 i ) + q(f * f - l) 

U) 

The resulting Euler -Lagrange equations are 


A + n • 0 

(5a) 

^ + G-^ p + (A * r) G^. - G^(A ' r)p 


+ (A • z) Vf G 2 = 0 

(5b) 

. T ~ 

0 i A • f = 0 

nr 

(5c) 

2r)J - ^ A = 0 

(5d) 


Since A, u, r, and v are all continuous except 
where sin intermediate boundary condition is im- 
posed (as will be shown later) , the three integrals 
are constant across staging points where continuity 
holds. However, for the oblate gravity model used 
in this analysis, only a single component of the 
above vector integral is constant, as can be 
verified by differentiation with respect to time. 

(A X v + JI X r) • z = constant (8) 

Transversality Equation 

The transversality equation for this problem 
is 

N _ t ± 

dJa V (C dt + A * dv + p * dr + a dm)^ 

1-1 

- <3m f (9) 

which is set equal to zero for an optimal solution. 
Reference 3 shows that: l) A and \i are continuous 

everywhere if there are no intermediate boundary 
conditions. If the intermediate boundary condition 
(assumed to occur at a staging point) is expressed 
as 

g(r,v) = 0 (10) 

reference 8_shows that the discontinuities in A 
and [i are ev^g and evpg, respectively. The vari- 
able e is used as an initial condition in the two 
point boundary value problem to satisfy the inter- 
mediate boundary condition (eq. (10)). 2) The 

equations that must be satisfied to optimize the 
duration of the powered and coast stages are de- 
rived in Ref. 3. The applicable results are pre- 
sented here. Let j be the first optimized powered 
stage. Then for constant jettison weight the equa- 
tion for optimizing stage t is 


x: (si - s? +1 ) = 0 (id 

i-j 


where o and f refer to initial and final values 
and the S functions are defined as 


S 


i 


C 

PT ’ 


a 


A * G + q ' 



6 




Pi J 0 

(l2a) 



(17) 


SjSO, P 1 => 0 (12b) 

where the right side of equation (12a) is obtained 
by using equation (7). For coasting stages 
=* 0) to be optimized, the equation 

C i =(A*G + JI-v) = 0 ( 13 ) 


must be satisfied for maximum payload. 

3) For free initial or final state variable x, 
the required or final condition for maximum pay- 
load (Ref. 4) is 



“ dr 

» ' 3£ 


(14) 


Initial Conditions 


If the Initial position and velocity are 
spec if I ed, _J;he initial values of any five of the 
six A and p may be used as variable initial 
conditions in order to satisfy the required final 
conditions of the two point boundary value pro- 
blem. In order to eliminate the difficulty asso- 
ciated with guessingjat values of the multipliers, 
the values of A and p can be expressed in terms of 
pitch and yaw attitude (\|r and 9) and rates U 
and 9). These equations majr be found in Ref. 4, 
appendix C. The values of A and p are then cal- 
culated from: 

% m XK , (l 5a ) 

^ - XK (15b) 

The value of A can be set equal to unity without 
loss of generality. The initial value of \ can 
be calculated in closed form, as will be shown 
by the following development . 

Final Conditions 

Final conditions for both the conventional 
and unconventional synchronous equatorial orbit 
mission require a circular orbit at synchronous 
orbit altitude with prescribed inclination. If 
the required inclination is non-zero, both the 
longitude of the ascending node and the injection 
point in the final orbit are free for optimiza- 
tion. As shown in Ref. 4, the corresponding 
auxiliary variational final conditions are 


(^gxv+^-gx?) • z - 0 

It will be shown later that equation (17) is 
satisfied for all functions g used herein. 


The calculation of A proceeds as follows : 
(A x v + p x r) * z = (AA x v) ■ z 
- (>£ x r) • z - (tf x 7) • z = 0 


A(A x v t A X r) * z 

A =■ jx =7 X 

(A x r) • z 


(18) 


Computing A with equation (18) guarantees that 
equation (l6a) will be satisfied. 


Intermediate Conditions 


As explained earlier, it is necessary to 
constrain the perigee radius at injection into 
the first parking orbit. Otherwise, the optimum 
solution would result in the parking orbit injec- 
tion and/or the equator crossing occurring at very 
low altitudes, thus violating spacecraft heating 
constraints. Therefore, the intermediate con- 
straint is 


g(r, v) = r p - r pjd = 0 

(19) 

where the desired value corresponds to the perigee 
altitude. Ey using equations found in Ref. 9, 
equation (19) can be written as 

rhr - r P> d = 0 

(20) 

where 


p s - (semi-latus rectum) 

G m 

(21a) 

e = 4/ 1 + (eccentricity) 

* G m 

* 

(21b) 

Q 

E = v ' v - JL (energy per unit mass) 

(21o) 

h = r x v (angular momentum per unit 

(21d) 


mass) 


A 

(x x V + n X r) • z = 0 

(16a) 

(Axv + ilxr) * (r x v) = 0 

( 16 b) 


If the desired inclination is zero, equations 
(16a) and (l6b) degenerate into one equation 
(zero inclination is equivalent to two final 
conditions, r * z = 0 and v • z = 0) , and only 
equation (l6a) must be satisfied. 

Since equation (l6a) Is a constant of the 
motion (eq. 8), it may be satisfied at the be- 
ginning of the trajectory, and used to calculate 
A . However, it must first be verified that jump 
discontinuities in A and p at intermediate boun- 
dary points do not change the value of the con- 
stant. This requires that 


The required gradients are calculated to be 


V- g 
v 


h(h x r - r p y ) 

» P 

e G m 

— (v x fi) - r 
* 


r 2 

-E 


(22a) 


(22b) 


It is easily shown that equation (17) is satisfied 
for the gradients in equation (22a) and (22b). 


In fact, equation (17^ is_satisfied for any 
function g of r, v, h and r * v. For such a 
function g, 


^ 6 




(h x r) 



r 


7 



V-g = ^r + 


and 


r ” dr 5h (7x ^ } (23) 


- dg - - 


^rgxv + TJpgxr^^vxv 


It should be recognized that there may be any 
number of fixed stages between tj and t k , etc. 

Also, the last three final conditions are evaluated 
at intermediate points in the trajectory. 

Appendix C 

Two Point Boundary Value Problem 


+ (h x r) x v + — r x v 
oh 3(r ■ v) 

+ |frxr+g|(vxh)x r 


+ y x r = 0 

a(p • v) 


(24) 


lienee the value_of Axv + ^ixris unaffected 
by the jump in A and \i resulting from a function 
g as defined above. 


Boundary Value Problem 

For the ATS-E mission, both fixed and 
optimum parking orbit coast times were considered. 
The transfer orbit coast time was always opti- 
mized, however, along with the durations of the 
first and second Centaur burns. Based on the 
preceding discussion of the transve reality equa- 
tion, the initial and final conditions for the 
two point boundary value problem are as follows 
for the case where the parking orbit coast time 
was optimized:. 

Initial Conditions 


t 

* ' 

9 

9 

e 

tj (First Centaur Burn) 
t^ (Parking Orbit Coast) 


t ^ (Second Centaur Burn) 
t m (Transfer Orbit Coast) 

(25) 


Final Conditions 


r^ ^ (Parking orbit) 
(r*z) = 0 
(v’z) = 0 

Z <<-=?*!> ■ 0 

I" j 


t-l - „ 

E <*< H > ■ o 

(A ’ G + |i * v) = 0 


If the desired final inclination is non-zero, 
then (r*z)_ and (v*z)_ = 0 are replaced by i, and 
(p x r + A x v) • (r x v) = 0. If the parking 
orbit coast time is fixed, then an Initial and 
final condition are removed. These are t^ and 



0 


The following technique was devised to syste- 
matically proceed from a simple, easily converged 
problem to the solution of the two point boundary 
value problem for a circular synchronous equatorial 
orbit. 

A trajectory is obtained to a slightly ellip- 
tical (parking) orbit with the desired perigee 
radius without plane change with a 90° launch azi- 
muth. This problem converges easily. Then the 
ascent burn time Is fixed at the value obtained 
and a variable length parking orbit coast, a fixed 
parking orbit perigee radius and a second burn are 
added. This problem Is targeted to the desired 
apogee and 180° argument of perigee for first 
equator crossing second burn. An inclination de- 
crease of about two degrees is then added to these 
final conditions and the problem is retargeted to 
the augmented final conditions. Now the transfer 
orbit coast (variable) and apogee burn (fixed or 
variable) are added. This trajectory Is integrated 
to the end with the converged initial guesses from 
the last step. The final conditions achieved will 
frequently be far from a circular synchronous equa- 
torial orbit. However, specify the final condi- 
tions actually achieved as the desired ones, and 
optimize the problem. The parking orbit coast, 
second burn, and transfer orbit coast durations 
will change. Now alter the achieved final condi- 
tions toward the desired ones judiciously in steps, 
retargeting at each step. In this manner, the 
desired final orbit conditions may be obtained. 

Now the ascent burn duration may be optimized. 

Any sizable change In a constraint or final condi- 
tion is best achieved by proceeding in steps. The 
problem is quite nonlinear. Attempts to plot 
initial conditions as functions of the final condi- 
tions for extrapolation purposes were made. They 
were generally unsuccessful. 
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Figure 1. - Conventional trajectory to circular synchronous 
equatorial orbit. 



Figure 2. - Circular synchronous equatorial orbit ascent profile. 


E-5732 



LAUNCH AZIMUTH, DEG 

Figure 3. - Separated spacecraft mass as a function of 
launch azimuth. 
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Figure 4. - Separated spacecraft mass as a function of 
final inclination. 
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Figure 5. - Transfer orbit inclination as a function of 
final inclination. 
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Figure 6. - Inertial velocity at apogee motor ignition as a 
function of final inclination. 
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Figure 7. - Latitude at second main engine start as a 
function of final inclination. 
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Figure 8. - Longitude at second main engine start as a 
function of final inclination. 
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Figure 9. - True anomaly at second main engine cut-off 
as a function of final inclination. 
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Figure 10. - Parking orbit coast time as a function of 
final inclination. 
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Figure 13. - Percentage of Centaur propellant burned 
during first burn as a function of final inclination. 
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Figure 14. - Final longitude as a function of final 
inclination. 
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